Binary black hole merger: symmetry and the spin expansion 
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We regard binary black hole (BBH) merger as a map from a simple initial state (two Kerr black 
holes, with dimensionless spins a and b) to a simple final state (a Kerr black hole with mass m, 
dimensionless spin s, and kick velocity k). By expanding this map around a = b = and applying 
symmetry constraints, we obtain a simple formalism that is remarkably successful at explaining 
existing BBH simulations. It also makes detailed predictions and suggests a more efficient way of 
mapping the parameter space of binary black hole merger. Since we rely on symmetry rather than 
dynamics, our expansion complements previous analytical techniques. 



In binary black hole (BBH) merger, two black holes 
A and B (with masses M a , Mb, and spins a, b) inspiral 
due to the emission of gravitational radiation and even- 
tually merge to form a final black hole with mass m, 
spin vector s, and recoil (or "kick" ) velocity k. How do 
the final quantities {m, k, s} depend on the initial quan- 
tities {M a , Mb, a,h}? This is a classic problem in gen- 
eral relativity (GR), with important implications for as- 
trophysics, cosmology and gravitational wave (GW) de- 
tection. For example, when two galaxies merge, their 
central supermassive black holes also merge. The final 
quantities {m, k, sjfrom these supermassive BBH merg- 
ers are linked (see [1| and references therein) to a variety 
of astrophysical observables including: (i) the quasar lu- 
minosity function; (ii) the location of quasars relative 
to their host galaxies; (iii) the orientation and shape of 
jets in active galactic nuclei; (iv) the correlation between 
black hole mass and velocity dispersion in the surround- 
ing stellar bulge; (v) the density profile in galactic cen- 
ters. The quantities {m, k, s} are also intimately related 
to the spectrum of quasi-normal ringdown modes after 
BBH merger — a key observable for probing black holes 
and strong-field GR with GW detectors. 

Following recent numerical breakthroughs 0, H, 3 , a 
number of groups can now simulate entire BBH merg- 
ers. In particular, they can choose a set of initial quanti- 
ties {M a , Mb, a, b} and compute the corresponding final 
quantities {m, k, s}. As more of these extremely time- 
intensive simulations have gradually accumulated, cer- 
tains patterns and trends have emerged. In previous 
work, some of these patterns have been described by em- 
pirical fitting functions which are loosely inspired by 
(but not derived from) post-Newtonian formulae Q. In 
this paper, we show how these same patterns (and others) 
may in fact be derived from elementary symmetry argu- 
ments. This perspective has several advantages which 
make it complementary to post-Newtonian and numeri- 
cal techniques. As we shall explain, the resulting formal- 
ism: (i) provides a simple conceptual understanding of 
BBH merger, accessible to non-experts; (ii) makes a host 
of new and derived predictions which go beyond the fit- 
ting formulae [5|, [(J ; (iii) suggests an efficient way to map 
out the parameter space of BBH mergers with simula- 
tions; and (iv) provides a map {M a ,Mb, a, b} — > {m, k, s} 



FIG. f: The orthonormal triad presented in the text. 



which is useful for astrophysical applications (includ- 
ing semi-analytic models or N-body simulations of black 
hole growth in galaxy mergers and dense stellar clusters) 
which wish to include BBH mergers, but cannot hope to 
follow the detailed merger dynamics. 

Although the merger process involves complicated non- 
linear dynamics, the initial and final states of the system 
are rather simple and symmetric. The initial state con- 
sists of two widely separated Kerr black holes, and the 
final state is a single Kerr black hole. The idea of this 
Letter is to see how much we can learn by considering 
the merger process as a map between these simple initial 
and final states, ignoring as much as possible the detailed 
dynamics in between. Under three operations (rotation 
"R," parity "P," and exchange "X"), the initial and final 
states transform according to simple rules. The depen- 
dence of the final state on the initial state is constrained 
to be consistent with these rules. We present a naive for- 
malism based on systematically applying these symmetry 
considerations to a well chosen Taylor expansion. 

Imagine two black holes, A and B, in a circular orbit 
fl7l |. The orbit gradually shrinks due to gravitational- 
wave emission, until A and B eventually merge. Con- 
sider an initial instant when the holes are far apart - 
far enough that they may be approximated as two Kerr 
black holes with well defined masses (M a and Mb) and 
dimensionless spins (a = S a /M% and b = Sb/Mfj), orbit- 
ing in a well defined plane that is perpendicular to the 
initial (dimensionless) orbital angular momentum Lq /M 2 



(where M = M a + M b , and G 



N 



= 1). Long after the 
merger is complete, the gravitational radiation has total 
energy -E ra d and total angular momentum J ra di an d the 
final Kerr black hole has dimensionless mass m = Mf /M, 



2 



dimensionless kick velocity k (relative to the center of 
mass), and dimensionless spin s = S//M?. 

At the initial instant mentioned above, we define an 
orthonormal triad {e^\e^ 2 \ e^ 3 )} as shown in Fig. [T] 
e^ 3 ) is the direction of the orbital angular momentum, 
is the direction from A to B, and e^ 2 ) = e^ 3 ) x e^ 1 ). 
Then the circular BBH's initial state is specified by 7 
numbers: the mass ratio q and the spin components 
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h = b • e (l > 



(1) 



relative to the orthonormal triad. Similarly, let us con- 
vert the final vectors into their separate triad components 



k • e« 



(2) 



If we apply a global 3-dimensional rotation R to the entire 
binary system (as if it were a single rigid body), initial 
and final quantities like a, b, k, and s rotate as vectors 
- as do the triad elements eW. Therefore, the corre- 
sponding triad components (a,, bi, k%, and s») transform 
as scalars under R. By working with triad components, 
all of our subsequent formulae are manifestly consistent 
with rigid 3-dimcnsional rotations R. 

We can view any final quantity / (such as m, Si, or ki) 
as a function of the initial quantities: 

/ = f(q,ai, a 2 ,a 3 ,bi,b 2 ,b 3 ). (3) 

Let us Taylor expand this function around a = b = 0: 

/ = f m ^ m2ms \ n ^ ns (q)a^ 1 a 2 n2 a^ 3 b^ 1 b^ 2 b'^ 3 . (4) 

Since a Kerr black hole has maximum spin Sj < Mf, |a| 
and |b| are both < 1 and it is not unreasonable to hope 
that the Taylor series might be convergent over most or 
even all of this range of initial spins. 

We will now use additional symmetries to restrict the 
coefficients in this expansion. First consider a parity 
transformation P that reflects every point of the binary 
system through the origin (the center of mass). Under 
P, the mass ratio q is unchanged, while the triad compo- 
nents transform as 



{ai,a 2 ,a 3 } -> {-oi, -a 2 ,a 3 } 
{bi,b 2 ,b 3 } -> {-61,-62,63} 



(5) 



Thus, each term on the right-hand-side of Eq. Q picks 
up a factor of (— l) 7 , where 



7 = m 1 + m 2 + n 1 + n 2 . 



(0) 



If / transforms under P as / — > (±) P /, the coefficients 
in Eq. ((4]) must satisfy the constraint 



/ 



m 1 m 2 m 3 \nin 2 n 3 



(q) = (±) P (-iP7 mi 



m 2 m 3 \n 1 n 2 n 3 
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(7) 

In other words, if / is even (odd) under P, then 
jm 1 m 2 m 3 |nir.2»3( g ) must van i sn w hen 7 is odd (even). 

Finally, we apply an "exchange transformation" X. 
This leaves the physical system absolutely unchanged, 
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TABLE I: Transformation under P, X, and PX, for various 
final quantities /. 



and simply swaps the labels of the two black holes, 
A <-> B. Under X, the mass ratio transforms as q — > 1/q, 
while the triad components transform as 



{ai,a 2 ,a 3 } 
{bi,b 2 ,b 3 } 



■bi, 
■ai, 



-b 2 ,b 3 } 
~a 2 , a 3 ] 



(8) 



If / transforms under X as / — > (±)^ /, the coefficients 
in Eq. (HJ must satisfy the constraint 

j"m 1 m 2 m 3 \n 1 n 2 n 3 _ (-|-) (_ \ yi f n l n 2n 3 \m 1 m 2 m 3 (X/g) 

(9) 

Equivalently, and more conveniently, if / transforms un- 
der PX (P followed by X, or vice versa) as / — > (±) PA - /, 
the coefficients in Eq. ([J| must satisfy the constraint 

ymim 2 m 3 |nin 2 n 3 _ ^-J-\ j^n 1 n 2 n 3 \m 1 m 2 m 3 (\ I q\ 

The transformation laws under P, X, and PX, for vari- 
ous final quantities /, are summarized in Table |TJ 

Eqs. (gl [TUTU]) imply a number of results that are exact 
(i.e. valid to all order in the Taylor expansion). For ex- 
ample, if the initial spin configuration has a oc b oc e^ 3 \ 
then 7 = and parity requires s oc e3 and k 3 = 0. In 
the equal-mass case, q = 1 = 1/q and Eq. (jTUJ) demands 
that all final quantities odd under PX (like fcj) vanish for 
a = b. Again in the equal-mass case, if the initial spin 
configuration satisfies (0,1,0,2, a 3 ) = (—b\,—b 2 ,b 3 ), then 
Eq. ijHJ) becomes the identity mapping and all quantities 
odd under X (like s\, s 2 ,k\,k 2 ) must vanish. 

Although these exact results are interesting, the real 
power of our formalism lies in the many approximate 
predictions that it makes. To illustrate this, consider 
the quantity k\. From Table [H k\ has (±) P = +1 and 
(±) PX = — 1. Start from the general Taylor expansion: 

fcj =fc™ 1, " 2m3|ni " 2 " 3 a™ 1 a™ 2 a™ 3 ^ 1 ^ 2 6™ 3 . (11) 



The zcroth-ordcr term k 
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vanishes by Eq. ([TO]). At 
first order, the Taylor expansion has six terms (one for 
each of the spin components {01,02,03,61,62,^3})- But 
constraints (J7J) and (fTU]) cut this down to a single term: 



, 0011000/ 

«i (a>3 



(12) 



Thus, at leading order, we predict k x oc (a 3 —b 3 ), a trend 
seen in simulations. But we also predict that this leading- 
order behavior should be corrected by quadratic terms of 
a specific form. Although the naive Taylor expansion 



FIG. 2: Equal-mass BBHs with a oc b cc e (3) . Final kick |k|, 
spin S3, and mass m are plotted versus a?,. The five curves in 
the top two panels correspond to the five simulation sequences 
in within each sequence 63 takes on the values shown in 
the bottom panel as 03 is varied. For presentation purposes 
we have exchanged 03 and 63 for the case 63 = —0.584. Cross 
and triangle data points for the mass m are taken from 
and [13], respectively. 

contains 21 new terms at 2nd order in spin, Eqs. |[7J) and 
pop reduce this to only 5 new terms, yielding the general 
2nd-ordcr formula 

, , 001|000/ , n . , 002|000/ 2 >2n 

k 1 = k 1 (03-63)+% (03-63) 

. , 200|000/ 2 .2\ , ,H0|000/ , , n 

+ k l (o 1 -6i)+fc 1 (0102-6162) 

. , 020|000/ 2 i2\ 1 ,100|010/ , , n Q \ 

+ & x (02-62) + % (a 1 6 2 -6 1 o 2 ). (13) 

This expansion could be continued to 3rd order and be- 
yond. At each order, we obtain more terms — although 
far fewer than a naive Taylor expansion would suggest. 

Our symmetry arguments are unable to determine the 
numerical values of the 6 coefficients (fc^ 01 ' 000 , etc.) in 
Eq. (fl~3"]) . These coefficients must be calibrated by 6 
equal-mass simulations with independent initial spin con- 
figurations. If each simulation measures a reliable value 
for fci, we can simply solve for these coefficients. Each of 
these simulations should be performed at a fixed value of 
some inspiral parameter (such as the orbital separation 
r$ or angular momentum Lq) which monotonically varies 
as the orbit shrinks. Having determined the coefficients, 
one can predict (with 2nd-order accuracy) the value of k\ 
resulting from any arbitrary configuration of the initial 
spins at the same fixed value of the inspiral parameter. 
Post-Newtonian methods 0] can then relate coefficients 
determined at other fixed inspiral parameter values 0. 

Similar arguments to those given in the k ± example 
apply to any final quantity /, including the final mass 



FIG. 3: Equal-mass BBHs with equal and opposite spins ly- 
ing in the orbital plane. Kick velocity £3, radiated angular 
momentum J ra d, final spin S3, and percentage of radiated en- 
ergy i? r ad are plotted against the angle cj> between the initial 
spin a and e' 1 '. The blue curves show fits to the square data 
points taken from while the red curves show fits to the 
triangle data points of flil j. 



m, and the magnitude and components of k and s. 

The explanatory power of this "spin expansion" for- 
malism is illustrated by a few simple examples. (For 
more details, we refer the reader to a subsequent paper 
Q, in which we discuss the formalism's new predictions 
in more depth, and test them in detail against currently 
available simulations.) We begin with the case of equal- 
mass BBHs with spins aligned or anti-aligncd with Lq 
{i.e. a cx b cx e^ 3 )). Expanding the final kick ki to 2nd- 
order in initial spins yields 3 terms for Ik} 2 , and 4 terms 
each for the final spin S3, and mass m [181 ]. Fig. [2] shows 
our best fits for this configuration. As seen in the figure, 
the data are well-described by the linear terms in the 
spin expansion, and the fits also show evidence for small 
second order corrections of the predicted form. 

Fig. [3] shows our fits for the "supcrkick" configuration 
(equal-mass BBHs with equal and opposite spins and 
a; = —hi = a(cos 0, sin 0, 0). The leading-order terms in 
the spin expansion explain the previously noticed [i~3l . fl4| 
behavior £3 = Aa cos((f) — (f>i) where A and (f>\ are con- 
stants (top panel in Fig. [3]). Keeping terms to 2nd- 
order, the spin expansion also correctly predicts that 
the final quantities / = { J ra d, S3, -E r ad} ah behave as 
/ = B + Ca 2 cos(20 — 2 ), where B, C, and are con- 
stants (bottom 3 panels in Fig. [3]). This cos(20 — <p2 ) 
behavior highlights the power of the spin expansion to 
go beyond previous linear post-Newtonian fitting formu- 
lae to uncover and explain new and essentially non-linear 
behavior in the simulations. 

Finally we examine the case of arbitrarily oriented ini- 
tial spins. From each of the 8 simulations in [l6| we have 
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FIG. 4: Equal-mass BBHs with equal spin magnitudes and 
generic spin orientations. Square data points 1 through 8 
indicate the final kick |k|, spin |s|, and mass m for eight sim- 
ulations listed from left to right in Table I of Red X's 
show our predictions for these final quantities. 



|s| and to. We use the coefficients obtained from our 
fits in Fig. [2] for the terms that depend only on 03 and 
63. Then our lst-order fits for |s| (with 3 free parame- 
ters) and to (with zero free parameters!) are shown in 
Fig. (H}. We stress that this data set is not described by 
any previous fitting formula. The spin expansion gives 
the first explanation for the distribution of points Fig. |4j 
In this Letter, we introduce a new spin expansion of 
final quantities / in triad components a, , 6; of the initial 
spins. Using the transformation properties of /, a,, bi 
and q under parity P and exchange X, we dramatically 
reduce the number of terms that one might naively ex- 
pect. Without resorting to the sophisticated machinery 
of numerical relativity and the post-Newtonian expan- 
sion, we obtain some detailed (and often new) quanti- 
tative understanding of the final state of BBH merger. 
This clarifies the separation between the non-linear dy- 
namics of Einstein's equations and our more elementary 



non-dynamical considerations. 

Our approach complements both post-Newtonian ap- 
proximations and numerical relativity Post-Newtonian 
methods provide accurate predictions in the weak-field 
inspiral regime, but break down during the later stages 
of the merger. By contrast, the symmetries under P and 
X implicitly hold through the entire merger. 

Only numerical relativity can model the late stages 
of the merger, but simulations remain computationally 
expensive. The spin expansion offers enormous compu- 
tational savings in mapping the 7-dimcnsional parameter 
space of BBH initial states {q, a^, bi}. Even 10 grid points 
along each direction would mean 10 7 simulations — a 
hopelessly large number. However, the values {m, k it s^} 
from 16 independent simulations determine the spin ex- 
pansion coefficients up to 2nd order at fixed mass ratio 
q. Then, with 10 grid points for q, a mere 160 simula- 
tions could map the space {q,a,i,bi}. Further reductions 
are possible if the q dependence of our coefficients can be 
identified analytically. The spin expansion may also help 
in identifying systematic errors with forbidden geometri- 
cal dependence on the initial spins. 

The spin expansion is useful for astrophysics and cos- 
mology — e.g. allowing BBH results from numerical rel- 
ativity to be efficiently included in simulations of cosmo- 
logical structure formation or black hole growth. These 
simulations cannot resolve the short scales relevant to su- 
permassivc BBH merger, and must instead rely on maps 
from the initial to final state such as those presented 
in this Letter. The spin expansion predicts these final 
quantities for arbitrary initial spin configurations, speci- 
fied once gravitational radiation (rather than dynamical 
friction) dominates the binary evolution. 
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